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Abstract. - We use the matrix product approach to construct all optimum ground states of 
general anisotropic spin-2 chains with nearest neighbour interactions and common symmetries. 
These states are exact ground states of the model and their properties depend on up to three 
parameters. We find three different antiferromagnetic Haldane phases, one weak antiferromag- 
netic and one weak ferromagnetic phase. The antiferromagnetic phases can be described as spin 
liquids with exponentially decaying correlation functions. The variety of phases found with the 
matrix product ansatz also gives insight into the behaviour of spin chains with arbitrary higher 
spins. 



Introduction. Haldane's conjecture Jl], || about the fundamental difference between 
integer and half-odd-integer spin chains has triggered a strong interest in low-dimensional 
systems with arbitrary spin S, both theoretically and experimentally. New materials have 
been found which can be considered as realizations of spin chains with S > i. For low- 
dimensional models mean-field theories are usually not very reliable due to the importance of 
quantum fluctuations. Therefore exact results become very important. They allow to explore 
possible phases and study their properties. Since only few models are solvable exactly, e.g. 
by Bethe- Ansatz ||, it is necessary to reach out for alternative methods. Next to numerical 
procedures like DMRG powerful analytical approaches are available. In this paper we 
use the idea of optimum ground states (|-||] which allows to construct systematically exact 
ground states. The method is not restricted to one-dimensional spin systems, but can also be 
generalized to construct ground states for quantum spin systems in arbitrary dimensions , 
Hubbard models or stochastic processes Q. 

A prominent example for a spin-2 chain is (2,2'-bipyridine)trichloromanganese(III) where 
the spin-2 is carried by Mn 3+ and the magnetic interaction by the Mn — CI — Mn orbital 
overlap. Experiments by Granroth et al. |0| show antiferromagnetic behaviour with expo- 
nentially decaying correlation functions. In the present letter we will construct states which 
exhibit exactly this behaviour with a non-degenerate ground state. Not only antiferromag- 
netic phases can be realized for a spin-2 chain but also weak antiferromagnetic and weak 
ferromagnetic ones. As the method of optimum ground states allows us to calculate ground 
state expectation values and properties for these states, it provides a good overview over the 
large variety of phases encountered in spin-2 chains and their properties. 

In the following we consider quantum spin chains described by a translational invariant 
Hamiltonian H = Y)u j\ h-ij with nearest neighbour interactions hij. 
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Definition. A global ground state \^o) is called optimum ground state (OGS) of H = 
^2(i j\ h-ij, if the global ground state energy Eq is just the sum of the local ground state energy 
values eo of the local Hamiltonians hij . 

In other words: If the local ground state conditions /iy|^o) = are complied (which can 
be achieved by adding the constant — eo to the local interaction hij) this yields the equivalence 

#|*o> =0 h ij \9 ) = 0, Vt,j (1) 

since then is a lower bound for Eq. The r.h.s. of (|I]) can be used as local conditions for 
the existence of an optimum ground state and its realization in terms of a matrix product 
ground state (MPG) (see below). Obviously OGS are states without finite-size corrections. 
In 0,ff3in the matrix product ground state approach has been described in more detail. 

For the spin-1 chain Q it is possible to construct one type of matrix product ground state 
for a Hamiltonian with the three symmetries 1) rotational invariance in the (x, y)-plane, 2) 
spin flip invariance, 3) translation and parity invariance. In the case of spin-2 chains with these 
symmetries and nearest neighbour interaction we find five different non-trivial MPGs. It can 
be understood easily that there exist no more than these five MPG except for exponentially 
degenerate ones Q. 

Usually it is most common to find exact or approximate ground states or even the full 
spectrum for a given Hamiltonian. The modus operandi we use here is the other way around. 
First we construct the optimum ground states for a given system and then we determine 
the subspace of Hamiltonians for which these are the ground states. Describing the local 
Hamiltonians hij in terms of local eigenstates \vk) of Sf + S*| (eigenvalue s) and the parity 
operator Pjj (eigenvalue p) and the spin flip invariance lead to — J2k ^k\vk){vk\ (see jig]). 
For some sets of quantum numbers s,p the \vk) are not determined completely by the imposed 
symmetries. Additional superposition parameters control the "orientation" of the basis states 
\vk) in their respective subspace. Requiring that a given OGS is an exact ground state of 
H leads to restrictions on the \k and the superposition parameters of the Iffe)- For a spin-2 
system with nearest neighbour interaction and above mentioned symmetries the total number 
of parameters is 22 pX( |. Two parameters, the energy off-set and a scale, are trivial which leads 
effectively to a 20 parameter model. For the most general model there are 25 local energy 
eigenstates. 

As expected such a model with 20 parameters in general has very complicated ground 
states of the global Hamiltonian H . The concept of optimum ground states and its realization 
through MPG yields a class of structural simple states. For these it is possible to calculate 
ground state expectation values for arbitrary operators. The isotropic valence-bond-solid 



(VBS) model |17| emerges as a special case and the properties of this model turn out to be 



generic for a more general subspace of the 20 parameter space. 

The matrix product ground state. Following the procedure explained in |(| one can 
define several matrices which are suited for constructing matrix product ground states of the 
spin-2 chain. Such a matrix consist of single-spin states at site i as their elements and the 
product of two such matrices is defined as a matrix product with the tensor product ® for 
the spin states 

( m W. m (^ = ^ m J) 8m ^), (2 ) 
k 
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Definition. A matrix product ground state (MPG) is a global ground state of a spin chain of 
length L, which for periodic boundary conditions can be written in the form 

|f )=tr(m< 1 ).m( 2 )-m< 3 )..-m< i )) , (3) 

where tr stands for "trace over the matrix space", i.e. (m*- 1 ) • ■ ■ ' TO ' L ') MAI - F° r a 
MPG to satisfy the above optimum ground state condition, hi.i+i\^o) = must hold for all 
i. Due to the product structure of the ground state |^o) t ms reduces to 

f m d) . m (*+iA = o for all i, (i, v. (4) 
This implies that each element of the product matrix 

( m (0 >TO (*+i)) is a local ground 
state (of hi^i + i). These conditions lead to several restrictions on the parameters from the 
most general model. So the matrix product ground state becomes a global ground state of the 
specified parameter subspace. The MPG itself turns out to depend on up to three parameters. 
To calculate expectation values and other physical properties a transfer matrix approach Jl4| 
can be used. Results thus obtained are presented in the following sections. 

MPGs on the spin-2 chain. - We have investigated systematically matrix product states 
on a spin-2 chain with unique or finitely degenerate ground states and periodic boundary 
conditions. We found five MPGs with non-trivial product structure and different properties. 
In the following only the most significant properties will be presented. A more complete 
account will be published elsewhere p5[ . 

Haldane-Antiferromagnet-A. - The first MPG is defined by the homogeneous product 
|^o)( a ) x t 7) = tr \ Ylf of matrices 



m w = 





0). yfr 


l). a 


Tr 






a 


2), V* 


T) 




(5) 



This representation with three continuous, real parameters a, a;, 7 uses the canonical spin-2 
basis states, i.e. the eigenstates of the S z -operator Sf\s z ). = s z \s z ). and S* ji 5 '). = — s z \s z ). 
at site i. Since all Zii^+i should be the same (for all i) we have to solve (Q) for one arbitrary 
i only. These conditions lead to a 12 parameter subspace and an additional trivial parameter 
for the scale. In general inequalities must hold for these parameters. For more details see fl5|| . 
This parameter space includes the isotropic point. The ground state itself depends on the three 
parameters a,x,j and is unique, which can be proven rigorously by complete induction ]l6| . 
For the parameter set x = —3, a = \/6,~f — — 2 (see fig. |l|) one obtains the isotropic point 
of the model where the matrix product ground state reproduces the VBS state E^]. The 
corresponding Hamiltonian can be written as sum of the projection operators onto the S z = 3 
and S z = 4 multiplets. For any other set of parameters the model has an anisotropy along 
the z-axis. However, in any case the ground state is unique as long as a ^ 0. For a — the 
matrix product Ansatz (MPA) is not unique; the degeneracy grows with system size as 3 L . 
This case is not considered in the following. 

The ground state defined by (||) is antiferromagnetic in the sense that all single-site mag- 
netisations vanish: (5 Z ) = (S x ) = (S v ) = 0. It follows that the square of the fluctuation is 
simply given by 

(AS 2 ) 2 = ((S z ) 2 ) - (s z ) 2 = ((S z ) 2 ) e [0,4]. (6) 
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Fig. 1 Fig. 2 

Fig. 1 - Longitudinal (thin) and transversal correlation length £ for x = —3 and 7 = —2. This special 
choice of parameter values includes the isotropic point for a = v6 where the correlation lengths are 
equal. 

Fig. 2 - For x = 0, 7 = 1 the correlation lengths diverge in the limit a — > 0. Note that the crossing 
point does not correspond to an isotropic Hamiltonian. 



The maximum (AS Z ) 2 = 4 is reached in the limit \a\ — > 00 and AS* 2 = for I7I — > 00. In the 
latter case the dominant contribution is proportional to 7 L |000- • -0). The 2-site correlation 
functions decay exponentially to zero. The longitudinal correlation function is given by 

(SfS?) = (SfSI) (sign (1 - a 2 )) r c-^/Z', (r > 2), (7) 

with the longitudinal correlation 

A = i ((1 + a 2 + 7 2 ) + y/(l + a 2 - 7 ) 2 + 8x 2 ) . (8) 

The first part of the function is the expectation value of the nearest neighbour correlation 
which is antiferromagnetic ((SfS^) < 0) and varies from (for 7 — > 00) to —4 (for a 2 — > 00) 
. For large a the correlation function alternates in r as an easy axis anisotropy is observed. 
The transversal correlation function decays exponentially, too, 

(SfS?) = (SfS?) (\x\ + 7 ) r e-(''- 2 >/«', = In 

Both correlation lengths and £t are finite (except for a = x = 0, see fig. ||) and so the 
model is non-critical. The a-dependence (see fig. [j] & ||) looks quite similar to the spin-1 MPG 
correlation lengths presented in |(| . For large a the anisotropy is of easy-axis type with > £ t 
whereas for small a it is of easy-plane type. Indeed this MPG is the analogue of the matrix 
product ground state of the spin-1 chain, e.g. the isotropic point (see fig. [j]) is included in 
this phase for both spin chains. For both chains the antiferromagnetic phases (including the 
two antiferromagnetic ones described in the next subsections) are Haldane phases with 1) a 
unique ground state, 2) exponentially decaying correlation functions, and 3) a gap to the first 
excited state. Here we only have shown 1) and 2) of this scenario, but the existence of a gap 
can be shown along the lines of the proof outlined in |l7]] . 

Haldane- Antiferromagnet-B. - Different from the spin-1 chain it is possible to construct 
further matrix product ground states for the spin-2 case. A simple one is a homogeneous 



in 



\x\ +7 



(9) 
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product of matrices 





o) 






( sfc 


T) 


a 


o) ) 



a e K,<r± 1, 



(10) 



which in contrast to the Haldane-antiferromagnet-A depends on one continuous parameter a 
and one discrete parameter a. The state |^o)(«) = tr (jjf mA looks quite similar to the one 
presented in || for the spin-1 case. Again it can be shown that it is a unique ground state 
of H in the appropriate subspace. Even the expectation values of all combinations of the Sf 
operators look the same as in the spin-1 case, only expectation values of S* 4 + or S~ operators 
vary by a factor. The properties of this phase are much the same of the one before and so we 
report some expectation values only. The single-site magnetisations vanish and the square of 
the fluctuation just as the longitudinal correlations function are the same as the ones for the 
spin-1 chain jq]. Different from the spin-1 chain the transversal correlation function reads: 



(S?S?) = -3\a\(sign(a)-a) 



C^lnll + H) 



(11) 



,i + M 

The correlation length £ t is the same as in H, only with a different prefactor. 

Haldane-Antiferromagnet-C. - Replacing in ( |l0|) the spin states |l) by |2) and |l) by 
1 2) we obtain a third matrix product ground state, which is the unique ground state in a 
16 parameter subspace of the most general 20 parameter model. The expectation values of 
S^-operators are the same as the ones above, except of a factor coming from the higher S z - 
eigenstates (| ± 2) instead of | ± l)). For other operators, e.g. S x = ^(S + + the method 

of transfer matrices shows that elements like (m^^\S + \mfi,o) become relevant. Because the 
matrix contains elements proportional to |2), |0) and |2) only, this always leads to zero and 
therefore ^SfS^) = 0. However, the expectation values of biquadratic operators of this type 
do not necessarily vanish, e.g. 

((SiSrf)=(^) 2 + 3(sign(a) + a)\a\( T ^ T - ] ) r . (12) 



1 



Weak- Antij err omagnet. - In the following we will show that also other antiferromagnetic 
phases can be realized by the matrix product technique, e.g. a phase characterized by a 
vanishing total magnetisation of two neighbouring sites (Sj + S^ +1 ) = 0, but with finite 
single-site magnetisation (Sj) ^ 0. If the sublattice is not fully polarized, i.e. \ (Sj)\ < 2, the 
corresponding state will be denoted as weak antiferromagnet. 

In order to realize a weak antiferromagnet using the MPA we introduce the two matrices 

(13) 

Using these matrices, two different MPGs can be constructed by assigning the matrices g and 
m to different sublattices: j^o^) = tr (ilf ^ m 2i-i • ff2i j and I^q 2 ^) = tr ^Jlf^ 2 92i-i ' m 2ij ■ 
In the first case to every odd lattice site a matrix m is attached and a matrix g to the even 
numbered sites. In the second case the situation is reversed. Therefore the ground state 
is twofold degenerate and each ground state depends on the parameter a. The single-site 
magnetisation alternates from lattice site to lattice site and can be written as 

\a\(x 2 

91 



(S* m ) = -{SI) = 1 + 



1) 



y/4 + a 2 (x 2 - l) 2 



(14) 
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For |a:| = 1 the sub-lattice magnetisation becomes (S^J = — \Sg) = 1. In the limit \x-a\ — > oo 
the magnetisation is (S z n ) = for values |x| < 1 and (S z n ) = 2 for \x\ > 1. For \xy/a\ — > oo 

but |a| -> a strict Neel order with |*^) = 1 2 22 . . . 2) and |*[, 2) ) = 1 222 . . . 2) is realized. 
The magnetisations in x- and y-direction vanish for the whole subspace 

(S x m ) = (SI) = (S* g ) = (S*) = 0. (15) 

The fluctuation AS Z varies in the range of to 1. It takes its maximum for |x| = 1 and 
\a\ — > oo. The longitudinal correlation function is given by 



(si%) = {-ir +i ((s* m ) 2 +A ie -^ 



) (r > 2) 



(16) 



and alternates in r. It decays exponentially to a constant contribution (S^) which is a 
result of translation invariance breaking. The amplitude Ai is positive and depends on the 
parameters a and x |l6| . 

The transversal correlation function decays exponentially to zero. Again, the correlation 
lengths stay finite for all values. In the limits a — ► 0+ and a — > 0— both correlation lengths, 
the transversal and the longitudinal, diverge. For a — the ground state degeneracy grows 
exponentially with the chain length. 

Weak-Ferromagnet. - The next two matrix product ground states look quite similar to 
the ones above. The main difference is that the translational invariance is not broken here. 
With the matrices 



1) 
0) 



2) 
1) 



the two MPG 



tr 



(n, 1 



and 




(17) 



can be constructed. Due to 



the spin flip symmetry 2) one has to ensure that both states are ground states of H which 
then is twofold degenerate. Just as expected the single-site magnetisation in z-direction reads 
(S z ) = —(S z ) = 1 and in x- resp. in ^-direction (S x ) — (S y ) = 0. The square of the 
fluctuation AS Z shows the same a-dependence as the one from the spin-1 MPG 



1 



(r > 2). 



The correlation lengths are the same as in the spin-1 case, 



l 



(18) 

(19) 
(20) 



Differences can be found in the longitudinal correlation function itself. The possibility to 



choose either j^o" 1 ') or I^q 3 ^) l ea ds to a finite contribution to the correlation function 



,(<?) 



(SfS z ) 



1 



(i-H) 2 Vi + H 



1- a 



(21) 



This partially polarized state is structurally similar to an antiferromagnetic state. It can be 
viewed as a Haldane-type state, but with finite magnetization. This is also reflected in the 
behaviour of the correlation function ( |2l| ) which agrees with that of the Haldane-AF-B up to 
the long-range order part which is zero in the latter case. 
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Conclusion. - The results presented here show that for a spin-2 chain the spectrum of 
possible MPG is much larger than for a spin-1 chain. We found three antiferromagnetic phases 
with unique ground state, exponentially decaying correlation functions and finite excitation 
gap. Therefore the corresponding phases can be classified as Haldane phases or spin liquids. 
Their structure is similar to that found in the case of the spin-1 chain Q. Similar antifer- 
romagnetic phases can not be constructed for the spin-| chain However, in this case a 
weak antiferromagnetic and a weak ferromagnetic phase exist. Related states can also be con- 
structed in the present spin-2 case. The weak antiferromagnet shows exponentially decaying 
correlation functions with long-range order. The ground state is twofold degenerate reflecting 
the breaking of translation invariance and leading to a finite sublattice magnetisation. The 
weak ferromagnet has a twofold degenerate ground state, but full translation invariance. The 
magnetisation per site in z-direction takes the constant value 1, which is just half of the fully 
polarised ferromagnet. Interestingly, despite the finite magnetisation, this state is structurally 
similar to an antiferromagnetic state. 

We believe that the exact results presented here are generic for spin-2 chains in a similar 
way as previous results (6[|l^] are for spin-1 chains. Here for the isotropic case, i.e. the 
bilinear-biquadratic chain, an exact solution using OGS is only possible for one point, the 
AKLT-chain (l7). However, a whole extended phase exists with the same properties. In this 
sense the exact results allow to study the generic properties of such phases without having to 
rely on approximate or numerical methods. 

Finally we like to point out that the results found here can be extended to even larger 
spins just as the Haldane-antiferromagnets-B and C can be viewed as analogues of the spin-1 
states. 
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